Swift monitoring of IGR J16418-4532 
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Abstract. We report on the Swift observations of the candidate supergiant fast X-ray transient 
(SFXT) IGR J16418-4532, which has an orbital period of -3.7 d. Our monitoring, for a total of 
—43 ks, spans over three orbits and represents the most intense and complete sampling along the 
orbital period of the light curve of this source. If one assumes a circular orbit, the X-ray emission 
from this source can be explained by accretion from a spherically symmetric clumpy wind from a 
blue supergiant, composed of clumps with different masses, ranging from —5 x 10 16 gto 10 21 g. 
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INTRODUCTION 

Supergiant Fast X-ray Transients (SFXTs) are a class of high-mass X-ray binaries with 
OB supergiant companions. Their hallmark is the occurrence of short (~hours) outbursts 
during which the luminosity can increase by 3-5 orders of magnitude (typically, up to 
10 37 erg s ). Either the clumpy structure of the wind from the companion or a centrifu- 
gal and/or magnetic gating have been suggested to be responsible for the outbursts. 

The hard-X-ray transient IGR J16418-4532 was discovered by INTEGRAL in 2003 
and later tentatively classified as an SFXT on the basis of its outbursting behaviour 
[1]. At soft X-rays XMM-Newton observations showed a heavily absorbed pulsar with 
Pspjj! ~ 1.25 ks [2]. An orbital periodicity of ~3.74 d was discovered with the Svvz/f/BAT 
and the RXTE/ASM [3, 4] and a hint of a total eclipse, consistent with a supergiant 
companion, has been observed in the BAT data. We analyzed all the Swift data collected 
on IGR J 164 18-4532, including a monitoring campaign carried out in 2011 July that 
spans over three orbital periods (see [5] for a complete log of the observations). 

RESULTS 

For the binary orbital period analysis we retrieved from the BAT Transient Monitor page 
[6] the 15-50keV BAT light curves (averaged over Swift's orbital period of ~5400s) 
covering the data range from 2005 February 12 to 2011 July 12. Fig. 1 (left) shows the 
epoch-folded BAT light curve for P orb = 3.73886(28) d and MJD 53560.20000 as r ep och 
[3, 4]. The BAT data show an eclipse centred at <j) ~ 0.55. The depth is consistent with a 
total eclipse, and the eclipse FWHM duration is A0 = 0.17 ±0.05 of the orbital period. 

Fig. 1 (right) shows the 0.3-10 keV XRT epoch-folded and background-subtracted 
light curve of the whole 2011 July campaign. In order to convert count rate to flux, 
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FIGURE 1. Left panel: BAT light curve folded at P orb = 3.73886 d and r epoch = MJD 53560.20000. 
Right Panel: XRT 0.3-10keV flux light curve of 1GR 116418^4-532, folded at the same period and epoch. 
The data were collected in 201 1 February and June (grey), and in 201 1 July (black). Downward-pointing 
arrows are 3(7 upper limits. The open circle is a detection obtained by combining three observations; it is 
the lowest point of the monitoring and corresponds to an unabsorbed 0.3-10keV flux of ^6.0 x 10~ 13 erg 
crrr 2 s -1 . The peak count rate is reached on MJD 55610.08, at 2.1 x 10~ 9 erg cirr 2 s~ ! . 



we adopted conversion factors derived from spectral fits of each observation (adopting a 
simple power-law model) and subsequent interpolation across observations where fitting 
was not feasible (see [5] for more details). The data start at phase 0.16 and cover three 
full periods. Superimposed on the long-term orbital modulation, which follows that seen 
with BAT, flaring is observed on short time scales (a behaviour typical of most SFXTs). 
In particular, the eclipse lasts A<j> ~ 0.2 and is centred at <j> ~ 0.55. 



DISCUSSION 

The Swift folded light curves of the source show a dip consistent with zero intensity. 
If interpreted as an eclipse, it can be used to infer the nature of the stellar companion. 
Excluding the eclipse, the BAT lightcurve does not show substantial flux modulation 
along the orbit, suggesting that the eccentricity is not large. This is in agreement with 
the observation that in many X-ray binary systems with nearly circular orbits and short 
orbital periods (< 4 days), the eccentric orbit produced by the explosion of the supernova 
is circularized on timescale of 10 6 yr (consistent with the evolution timescale of an O 
star), due to viscous tidal interactions. Thus, for simplicity, we assume a circular orbit 
Cforb = 3.73886 d). Let us adopt the typical values of stellar masses and radii [7] for a 
primary of spectral type 08.5 and luminosity classes I, III, or V, and for the neutron 
star a mass of Mns = 1.4M , and a radius of 7?ns = 10km. By combining Kepler's 
third law with the relationship given by Rappaport and Joss [8] for the radius of the 
primary R, (R/a) 2 = cos 2 / + sin 2 / sin 2 e , where a is the semi-major axis of the orbit, i 
is the inclination, and 6 S is the eclipse semi-angle, we find that the observed A(j) = 0.2 
is inconsistent with a luminosity class III or V primary. We conclude that indeed the 
primary is an 08.51 supergiant star (with mass, radius, luminosity, effective temperature 
given by Martins et al. [7], M b = 31.5 M , R 0B = 21.4 R , logL B/£© = 5.65, 
r eff = 32274 K) and that the distance to IGR J16418-4532 is -13 kpc [9]. 



The XRT light curve shows a large dynamic range, up to 1-2 orders of magnitude 
within a single snapshot and of at least 3 orders of magnitude overall (even when the 
low-flow flux measured during the eclipse is excluded). While this dynamic range is 
lower than that observed in the classical SFXTs, which reach 4-5 orders of magnitude 
in dynamic range, it is indeed typical of intermediate SFXTs (e.g. IGR J18483— 031 1). 
The 201 1 February flare (at phase ~0.27 in Fig. 1), in particular, reached ~2x 10~ 9 erg 
cm~ 2 s~ ! (2-lOkeV, unabsorbed), corresponding to a luminosity of ~ 4 x 10 37 erg s _1 
at 13 kpc, which is typical of SFXT outbursts. 
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FIGURE 3. Observed and computed cu- TABLE 1. Wind parameter values for 
mulative luminosity distribution; see text. IGR J16418-4532. 

To reproduce the observed X-ray variability, we applied a clumpy wind model [10], 
which we briefly summarize here. This model assumes that the wind of OB supergiants 
is inhomogeneous, composed of dense clumps surrounded by a hotter and homogeneous 
wind. Clumps follow a power-law mass distribution, p(M c \) = fc(M c i/M a ) _ ^, where 

k = (1 — QM a /(M h — M a s ) is the normalization constant, M c \ is the mass of a 
clump, and M a and Mb define the extrema of the mass range. Assuming a spherical 
geometry for the clumps, the initial clump dimension distribution is given by Nm c1 
RL clumps s _1 , where R c \ is the radius of the clump. The clump mass loss rate M c \ is 
related to the total mass loss rate M tot by means of / = M c i/M tot , the fraction of wind 
mass in the form of clumps. Clumps and the homogeneous wind follow the ^-velocity 
law v(r) = Vco (1 -0.9983 j R O bA) /3 [11], where /3 is a constant in the range ~0.5 — 1.5, 
7?ob is the radius of the supergiant and Voo is the terminal wind velocity. 

Each Swift snapshot has a duration of ~ 1 ks, for a total of up to ~ 5 ks per day, which 
is the typical timescale for the flare durations observed in SFXTs. For this reason, we 
cannot establish for IGR J 164 18-4532 the number of observed flares and their duration, 
therefore we cannot use this information to compare the observed lightcurve with that 
calculated with the clumpy wind model. However, it is possible to compare the observed 
cumulative luminosity distribution out of eclipse with that calculated with the model. 
For the former, the SwiftlXKY count rate, measured in the 0.3-10 keV, was converted to 
the 0. 1-100 keV luminosity using the average spectral parameters obtained by [2, 12], 
and for a distance of 13 kpc. For the latter, we adopted the physical parameters for the 
08.51 primary and the neutron star secondary reported above, and reasonable the force 
multiplier parameters for a 08.51 star [see 5]: k = 0.375, a = 0.522, 8 = 0.099. 



Fig. 3 shows the XRT cumulative luminosity distribution (solid line) with its 90 % un- 
certainties (dashed lines) and also the cumulative luminosity distribution (red dot-dashed 
line) calculated with the clumpy wind model by adopting the following parameter val- 
ues: mass loss rate A/ to t = 6 x 10~ 7 M yr -1 , terminal velocity Voo = 1000 km s , j8 = 1, 
fraction of mass lost in clumps / = 0.9, mass distribution power-law index £ = 1.1, 
power-law index of the initial clump dimension distribution 7 = —1, minimum clump 
mass M a = 5x 10 18 gand maximum clump mass Mb = 10 20 g. Further acceptable solu- 
tions can be found with the wind parameters in the ranges reported in Table 1. 

The wind parameters obtained for IGR J16418— 4532 are roughly in agreement with 
those typical of HMXBs, with the exception of / and Voo. The obtained low value of the 
terminal velocity can be explained as follows: when the wind is highly ionized by the 
X-ray photons emitted by the neutron star, the high ionization modifies the dynamics 
of the line-driven stellar wind of OB supergiants (especially in close binary systems), 
leading to a reduction of the wind velocity in the direction of the neutron star [e.g. 13]. 
The mass loss rate is lower than the typical mass loss rate of 08.51 stars, which is of the 
order of m 2 x 1O~ 6 M yr _1 [14], but in agreement with the hypothesis that the mass 
loss rates derived from homogeneous-wind model measurements with the Ha method 
are overestimated by a factor of 2-10 if the wind is clumpy [e.g. 1 5] . However, we cannot 
exclude the possibility that the X-ray variability of IGR J16418— 4532 is not totally due 
to the accretion of an inhomogeneous wind. In fact, other mechanisms could be at work, 
like e.g. transient accretion discs and intermittent accretion flow onto a neutron star. 
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